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ABSTRACT 



Context. X-ray imaging observatories have revealed hydrodynamic structures with linear scales ~ 10 kpc in clusters of galaxies, such 
as shock waves in the 1E0657-56 and A520 galaxy clusters and the hot plasma bubble in the MKW 3s cluster. The future X-ray 
observatory IXO will resolve for the first time the metal distribution in galaxy clusters at the these scales. 

Aims. Heating of plasmas by shocks and AGN activities can result in non-equilibrium ionization states of metal ions. We study the 
effect of the non-equilibrium ionization at linear scales < 50 kpc in galaxy clusters. 

Methods. A condition for non-equilibrium ionization is derived by comparing the ionization time-scale with the age of hydrodynamic 
structures. Modeling of non-equilibrium ionization when the plasma temperature suddenly change is performed. An analysis of 
relaxation processes of the FeXXV and FeXXVI ions by means of eigenvectors of the transition matrix is given. 
Results. We conclude that the non-equilibrium ionization of iron can occur in galaxy clusters if the baryonic overdensity 5 is smaller 
than 1 1.0/t, where t <k 1 is the ratio of the hydrodynamic structure age to the Hubble time. Our modeling indicates that the emissivity 
in the helium-like emission lines of iron increases as a result of deviation from the ionization equilibrium. A slow process of helium- 
like ionic fraction relaxation was analyzed. A new way to determine a shock velocity is proposed. 

Key words. Galaxies: clusters: general; Atomic processes; Shock waves 



1. Introduction 



Clusters of galaxies are gravitationally bound structures of mass 
~ 10 14 - 10 15 Mq and size ~ 1 - 3 Mpc (for a review, see 
IKaastra et al.l 2008). Their mass budget consists of dark matter 
(~ 80%), hot diffuse intracluster plasma (< 20%) and a small 
fraction of other components such as stars and dust. The mean 
baryonic overdensity in galaxy clusters equal to 5 — nu/nu — 
200, where «h and «h are the mean hydrogen densities in galaxy 
clusters and in the Universe, respectively. 

Many chemical elements reside in galaxy clusters. The 
plasma temperatures kT^ 3-10 keV in galaxy clusters are close 
to the values of the K-shell ionization potentials of heavy ele- 
ments (Iz = Z 2 Ry, where Z is the atomic number and Ry the 
Rydberg constant). Emission lines from heavy elements were 
detected by X-ray telescopes from galaxy clusters. The cur- 
rent instruments (XMM-Newton, Chandra and Suzaku) have 
largely enhanced our knowledge on the chemical abundances of 
many elements. Metal abundances around 0.3 in Solar Units of 
lAnders & Grevesse (1989) were derived under the assumptions 
of collisional ionization equilibrium (for a review, IWerner et al.l 
2008). 

Non-equilibrium processes such as non-equilibrium ioniza- 
tion and relaxation of the ion and electron temperatures are usu- 
ally taken into account only in the outskirts of galaxy clusters 
and in the warm hot intergalactic medium (WHIM) where bary- 
onic overdensity 6 is less than 200 (e.g. Yoshikaw a & Sasakil 
2006; Pr okhorovl 2008). However, we show that non-equilibrium 
ionization can also be produced as the result of merging pro- 
cesses and AGN activity in galaxy clusters where the baryonic 
overdensity 5 > 1000. 
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Evidences for merging processes of galaxy clusters and AGN 
activity, such as shocks and hot plasma bubbles, were revealed 
by means of Chandra high-resolution observations. For example, 
strong shocks in the 1E0657-56 and A520 galaxy clusters propa- 
gating with a velocity of 4700 km/s and 2300 km/s, respectively, 
were derived by Marke vltch et al.l (2002) and Marke vitch et al.l 
(2005). The corresponding Mach numbers of the shocks are 3.0 
and 2.1. Hot plasmas inside bubbles arising from AGN activ- 
ity were detected in galaxy clusters (e.g. IMazzotta et aTl 2002). 
Heating of plasmas produced by shocks or AGNs can result in a 
non-equilibrium ionization state. 

In this paper we study the effect of non-equilibrium ioniza- 
tion near merger shock fronts and in hot plasma bubbles. We give 
a theoretical analysis of collisional non-equilibrium ionization 
in Sect. 2. We show numerically the importance of this effect in 
galaxy clusters in Sect. 3. We analyze helium-like and hydrogen- 
like non-equilibrium ionization states by means of eigenvectors 
of the transition matrix in Sect. 4. A new approach to determine 
the value of the shock velocity is considered in Sect. 5 and our 
results are discussed in Sect. 6. 



2. A condition for non-equilibrium ionization. 

Non-equilibrium ionization is often assumed in supernova rem- 
nants (e.g. IGronenschild & Mewel 1982. iMasail 199 4) and may 
be important in the WHIM dYoshikawa & Sasakil 2006). The 
non-equilibrium ionization state in the linked region between 
the Abell 399 and Abell 401 clusters was also studied by 
lAkahori & Yoshikawal (2008). We are going to show that non- 
equilibrium ionization can occur not only in the outskirts of 
galaxy clusters but also in galaxy clusters, in which merging pro- 
cesses and AGN activity play a role. We derive here a condition 
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on the baryonic overdensity for deviation from collisional ion- 
ization equilibrium, by comparing the ionization time-scale with 
the age of hydrodynamic structures. 

The number of collisions between electrons and an ion 
resulting in electron impact ionization per unit time is v = 
0"ionVthr", where cn on is the characteristic value of the ionization 
cross-section, Vthr corresponds to the threshold velocity (energy) 
of the ionization process, and h is the number density of elec- 
trons which have sufficient energy for electron impact ionization. 
The ionization time-scale is given by t lon = 1/v, and therefore 

1 



fion = -• (1) 

0"ionV t hr« 

As it was noted by lYoshikawa & Sasakil (2006) and 
lAkahori & Yoshikawa (2008), helium-like and hydrogen-like 
ions are interesting for the analysis of non-equilibrium ioniza- 
tion. Therefore, electron impact ionization of a helium-like ion 
will be considered in this section as a physically important case 
(a consideration of hydrogen-like ions is analogous). 

For helium-like ions the characteristic value of the ionization 



cross-section is approximately (see Bazylev & Chibisov 1981) 



(2) 



where ao = ft 2 /m e e 2 is the Bohr radius, Z is the atomic number. 

Since the ionization potential of a He-like ion is approx- 
imately Iz ~ Z 2 m e e 4 /T 2 /20, the electron threshold velocity 
Vthr = V2/z/'«e can be approximated as 



Vthr 



e 2 
Z— . 

h 
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The number density of electrons with energies higher than 
the ionization potential of a He-like ion is 
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where «o is the plasma number density, = yfZIzJkT is 
the dimensionless threshold momentum and fu{x) = yfT/n x 
exp(-x 2 /2) is the Maxwellian distribution. 

If the dimensionless threshold momentum p± r > 1 then a 
simplified form of Eq. is given by 



5 Vi n ° eXP l kT) 
Using Eqs. ©, ([3]) and (f5]) we rewrite Eq. ([D as 
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It is most convenient to write the plasma number density in 
terms of the baryonic overdensity «o = SClbPmt/mp, where the 
critical density is p cl ; t = l>H 2 l(%nG), and to denote the ratio of 
the thermal energy kT and the ionization potential by A — kT /Iz- 
Thus, 



(7) 



2ViZ m e V G ®(A) 
t,on * J^-VH 2 ^ S~ 

1 The exact value of the ionized potential of a He-like ions is obtained 
by changing Z with Z e ff = Z - 5/16. We consider the case Z » 1 and, 
therefore, the approximative value is sufficient. 



where <E>(T) = VI exp (1 /A). 

Ionization states will be non-equilibrium if the ionization 
time-scale t lon is longer than the hydrodynamic structure age 
t = tH~ 1 , i.e. f; on > tH~ x . This condition is equivalent to the 
inequality 



5< 
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In an important case of the iron ions (Z=26), the numerical 
value of the first dimensionless term on the right-hand side of 
Eq. (Ull is 



2 yj7rZ 3 mlm p e 2 G 
3 V2>i 3 // fi b 
and, therefore, 
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In rich galaxy clusters with plasma temperature of kT p \ ss 5 
keV, the value of the function <&(kT p i/Iz) is * 4.6. Therefore, in 
this case we find 6 < 11.0/t and conclude that, if the hydrody- 
namic structure age is of order 10 7 - 10 8 years (i.e. t lies in the 
range 6.6 x 10~ 4 < r < 6.6 x 10 3 ) then non-equilibrium ion- 
ization occurs in galaxy clusters where the baryonic overdensity 
6* 1000 (seeEq.[p . 

IMazzotta et al. 1 (2002) have estimated the age of a hot plasma 
bubble of diameter ~ 50 kpc to be t as 2.7 x 10 7 yr, which is much 
shorter than the age of the MKW 3s cluster. In the 1E0657-56 
and A520 clusters the downstream velocities of the shocked gas 
flowing away from the shock are 1600 km s~' and 1000 km s _I 
dMarkevitch et al.l 2002, IMarkevitch et all 2005), therefore the 
shocked gas covers a distance 50 kpc in times 3.1 x 10 7 and 
4.7 x 10 7 yrs respectively. Thus, in light of the above conclusion 
non-equilibrium ionization can occur at linear scales < 50 kpc 
in galaxy clusters in which merging processes and AGN activity 
present. 



3. Modeling of non-equilibrium ionization. 

Non-equilibrium ionization occurs when the physical condi- 
tions of the plasma, such as the temperature, suddenly change. 
Shocks, for example, can lead to an almost instantaneous rise 
in temperature and to a deviation from ionization equilibrium. 
However, it takes some time for the plasma to respond to an 
instantaneous temperature change, as the ionization balance is 
recovered by collisions. 

In this section we consider the following situation: the 
plasma temperature instantaneously increases from kTi =3.4 
keV to kT2 = 10.0 keV. Such a temperature change may cor- 
respond to a temperature jump at a shock with a Mach number 
M = 2.6 or to plasma heating by AGN activity. We assume that 
the age of the hot plasma region is 3 x 10 7 yr and the baryonic 
overdensity is 5 = 4000 which corresponds to the plasma num- 
ber density in the post-shock region in the A520 cluster (see Fig. 
2b of Markevitch et al. 2005). Following |Markevith| (2006) we 
assume that the electron and ion temperatures are equal. 

At the temperature kT\ =3.4 keV the ionic fractions of Li- 
like, He-like and H-like ions of Fe relative to the total Fe abun- 
dance are » 12%, * 77%, and * 9% respectively. Therefore we 
consider below four iron-ion-states Fe(+23), Fe(+24), Fe(+25) 
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Fig. 1. Dependence of the He-like (solid line) and H-like (dashed 
line) ionic fractions of iron on the dimensionless time f/(3 x 10 7 
yr). 



and Fe(+26). In the case the collisional ionization rate equation 
for each element is written as 



riQdt 
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(11) 



where n is the vector with four components (n+23, n+24, «+25, 
n+26), normalized such that 2?=23 «+i = 1, which correspond 
to the four iron-ion-states mentioned above, 7 +z and R +z rep- 
resent the rate coefficients for ionization and recombination 
from an ion of charge z to charges z+1 and z-1, respectively. 
All the coefficients necessary to calculate the direct ionization 
cross sections are taken from Arnaud & Rothenflug (1985), the 
radiative recombination rates are taken from IVerner & Ferlandl 
(1996), and the dielectronic recombination rates are taken from 
IMazzotta et al. 1 (1998). To solve the system of equations (fTTT i we 
use the fourth order Runge-Kutta method. 

The time-dependence of the He-like and H-like ionic frac- 
tions of iron is shown in Fig.Q] 

At the temperature kT2 = 10.0 keV the equilibrium ionic 
fractions of He-like and H-like ions of iron are « 27% and 
ss 45% respectively. Therefore, the He-like ionic fraction which 
equals 35% at time 3 x 10 7 yr does not reach its equilibrium 
value and non-equilibrium ionization occurs. However, the H- 
like ionic fraction almost reaches its equilibrium value at time 
3 x 10 7 yr. 

We now show that the effect of non-equilibrium ionization 
on the helium-like emission lines of iron can be significant and 
that non-equilibrium ionization leads to the increase of volume 
emissivity in the helium-like spectral lines. 

The helium-like volume emissivity for a chemical element 
of atomic number Z is given by 



fZ = "e"HA Z X («+(Z-2)G+(Z-2) + «+(Z-l)«+(Z-l)) 



(12) 



where « e is the electron number density, «h is the hydrogen 
number density, Az is the abundance of the considered chemi- 
cal element, «+(z-2) and «+(z-i) are the ionic fractions of helium- 
like and hydrogen-like ions respectively, Q+(z-2) is the impact 
excitation rate coefficient and ar+(z-i) is the rate coefficient for 
the contribution from radiative recombination to the spectral 
lines. Excitation rate coefficients are taken from lProkhorov et al.l 



(2009). Let us note the reduced volume emissivity in the iron 
helium-like emission lines as 



U 



£Z=26 

n e n H A z T 



(13) 



where F = Z~ 4 na^ yjlzhn e corresponds to the characteristic rate 
coefficient value (see also |Prokhorov| 2009). 

In Fig.[2]the reduced emissivity U when the ionic fractions 
are in ionization equilibrium is shown in the range of tempera- 
tures between 3.5 keV and 1 1 keV. 
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Fig. 2. Dependence of the equilibrium reduced iron volume 
emissivity in the helium-like lines on the plasma temperature. 



In the situation considered above we find that at temperatures 
kT\ = 3.4 keV and kTi = 10.0 keV the equilibrium values of the 
reduced volume emissivities are approximately equal. However, 
in the presence of non-equilibrium processes approximate equal- 
ity of these volume emissivities does not hold. Since the frac- 
tion of electrons with energy higher than the impact excitation 
threshold E ex m 6.7 keV is 27% at the temperature kT\ = 3.4 
keV and is much less than 72% that is at temperature kTi = 10.0 
keV, more effective impact excitation should be at temperature 
kT2 = 10.0 keV. Furthermore the non-equilibrium ionic frac- 
tion of helium-like iron in the region of temperature KT2 = 10.0 
keV is higher than the equilibrium ionic fraction (see FigQ]) and, 
therefore, non-equilibrium ionization leads to the increase of 
volume emissivity in the helium-like spectral lines. 

Using the dependence of the ionic fractions of iron on the 
dimensionless time f/(3 x 10 7 yr) (see Eq.[TT1i we study the time 
evolution of the reduced volume emissivity in the iron helium- 
like emission lines. This time evolution is shown in Fig. [3] 

Since the maximal value t/ max of the equilibrium reduced 
emissivity U(T) is at the temperature kT ss 6 keV and equal to 
=s 1.25 (see Fig. O, we conclude that the non-equilibrium value 
of U(t) (see Fig. [3]) is higher than the maximal equilibrium value 
U mdx up to the time t — 2 x 10 7 yr (up to the dimensionless time 
equal to 0.65). The value U eq of the equilibrium reduced emis- 
sivity at the temperature kT = 10.0 keV equal to * 0.8 is less 
than the non-equilibrium value of U (f) during the time interval 
t — 3 x 10 7 yr. Therefore, the iron abundance Az, eq derived from 
the assumption that ionization states are in equilibrium will be 
higher than the correct iron abundance value Az (see Eq. [T3l . 
which is given by 



A 7 =A 



Z,eq 



u eq 



(14) 
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Fig. 3. Evolution of the reduced iron volume emissivity in the 
iron helium-like lines in the region with temperature 10.0 keV. 
The dimensionless time is given by t/(3 x 10 7 yr). 



The eigenvector which corresponds to the eigenvalue Aq de- 
termines equilibrium ionic fractions at temperature kT = 10.0 
keV. 

Since A3 is the smallest absolute value of the eigenvalues (ex- 
cluding /lo which does not correspond to any relaxation process) 
the process which corresponds to the eigenvector V3 is the slow- 
est (see Eq. UTT l. This slow process corresponds to the increase 
in the FeXXVII ionic fraction due to decreases in the FeXXV 
and FeXXVI ionic fractions. However, the absolute value of the 
second component of V3, which corresponds to the decrease in 
the FeXXV ionic fraction and equals 0.61, is higher than the ab- 
solute value of the third component of V3, which corresponds 
to the decrease in the FeXXVI ionic fraction and equals 0.14. 
Therefore, the variation in the helium-like FeXXV ionic frac- 
tion which is proportional to the value of the second component 
of V3 is more substantial during this relaxation process than the 
variation in the hydrogen-like FeXXVI ionic fraction. 



4. An analysis of He-like and H-like non-equilibrium 
ionization states by means of eigenvectors 

In the previous section we showed that the ionic fraction of He- 
like iron ions can remain in non-equilibrium while the H-like 
iron ionic fraction almost achieves equilibrium. This somewhat 
paradoxical behavior can be more easily understood by means 
of eigenvectors of the transition matrix M, which is (see Eq. fTTb 



M = « f ag e 
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(15) 
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where f age = 3 x 10 7 yr is the age of the hydrodynamical structure 
(see Sect. 3). 

Here we calculate the values of the eigenvalues of the tran- 
sition matrix and the corresponding eigenvectors, and show 
how the ionic fraction of He-like iron ions can remain non- 
equilibrium longer than that of H-like iron ions. 

The eigenvalues A of the transition matrix M are derived 
from the equation 



Det(M - A X E) = 0, 



(16) 



where E is the unit matrix. 

One of the eigenvalues of the transition matrix M is of the 
form Eq. dl~5b equal to zero (Ao = 0). Consequently, the ioniza- 
tion equilibrium is achieved in the end. 

The solution of the system of differential equations (Eq. fTTT i 
can be written as 



n = V CjViexpl^-l, 

i=o 



(17) 



where c\ are constants, n is the vector («+23, n+24, «+25, «+2e) and 
Vi are the eigenvectors of the transition matrix M. 

At the temperature kT = 10.0 keV we derive three eigenval- 
ues which equal to A\ ~ -17.80, A 2 ~ -3.40 and A3 ~ -1.46. 
The eigenvectors which correspond to the derived eigenvalues 
are respectively 
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Fig. 4. Dependence of the He-like (solid line), H-like (dashed 
line) and fully ionized (dot-dashed line) ionic fractions of iron 
on the dimensionless time f/(3 x 10 7 yr). 



The quantitative difference between the variations in the He- 
like and H-like ionic fractions is illustrated in Fig. [4] 



5. Mach number determination 

Clusters of galaxies form via gravitational infall and mergers of 
smaller mass concentrations. In the course of a merger, a signif- 
icant portion of the kinetic energy of the colliding subclusters, 
that carried by the gas, is dissipated by shocks. 

The shocks in the A520 and 1E0657-06 clusters have Mach 
numbers M = 2-3, derived from the Rankine-Hugoniot jump 
conditions (for a review, see IMarkevitch & Vikhlinin 2007), re- 
lating the density and temperature jumps at the shock and the 
Mach number, M = v/ci, where c\ is the velocity of sound 
in the pre-shocked gas and v is the velocity of the pre-shock 
gas in the reference frame of the shock. Thus, if the pre-shock 
and post-shock temperatures (7"i and To, respectively) are deter- 
mined from observations, the Mach number of the shock can be 
derived from the equation: 



II 
Ti 



(2yM 2 sh -(y-l))x((y-l)Ml + 2) 

(y+l) 2 /< 



(19) 



where y is the adiabatic index. It is usually assumed that the pre- 
shock velocity in the reference frame of the shock is equal to the 
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shock velocity in the reference frame of the galaxy cluster and 
that the adiabatic index is y — 5/3 (see Markevit ch & Vikhlininl 
2007). 

The comparison of the X-ray image and gravitational lensing 
mass map of the 1E0657-06 merging cluster (IClowe et al. 1 2006) 
shows that the mass peak of the subcluster is offset from the 
baryonic mass peak. IClowe et al. 1 (2006) interpret this as the first 
direct evidence for the existence of dark matter. 

Such merging clusters offer the unique opportunity to study 
gas physics through direct comparison of the observed shock 
properties with the predictions of gas + dark matter modeling 
(e.g.lProkhorov & Durretl2007;|Springel & Farrar|2007). In this 



section we provide a new way to derive shock parameters based 
on measurements of the flux ratio of the FeXXV and FeXXVI 
iron lines. 

The fluxes of the FeXXV and FeXXVI lines have the same 
dependence on the metal abundance, as well as on the emission 
measure, their ratio is independent of these parameters. This iron 
line ratio can therefore be used to determine the temperature of 
the intracluster gas (e.g. Nevalaine n et al.1 2003) and the presence 
of supra-thermal electrons (e. g. IProkhorov et aLl 2009). 

Taking into account both electron-impact-excitation and ra- 
diative recombination the iron line flux ratio is given by 



R 



»+246 p e xxv + n +25ffR R 2 Fe XXV 
"+25 Q 



FeXXVI 



+ n 



1-2 ' 
-26<* RR FeXXVI 



(20) 



where the rate coefficients are Q 



1-2 



2 Z ^ FeXX\ 

a b(<a) 

Gftxxvi = 22 Spexxvi^- and S F * X ' XVI are the impact- 

a b(<a) 

excitation rates. The excited states b correspond to the upper 
levels of the He-like triplet and the H-like doublet, and the ra- 
diative, branching ratios are given by, 
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Fig. 5. Evolution of the iron line flux ratio in the region with 
temperature 10.0 keV. The dimensionless time is given by t/(3 x 
10 7 yr). 
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and a RR FeXXV and a RR FeXXVI are the rate coefficients for the 
contribution from radiative recombination to the spectral lines 
FeXXV (He-like triplet) and FeXXVI (H-like doublet), respec- 
tively, and Aac are the transition probabilities. 

Below we study the situation which was considered in Sects. 
3 and 4. The variation of the iron line flux ratio in the region 
with temperature 10.0 keV as a function of dimensionless time 
is shown in Fig. [5] 

Let the downstream velocity of the shocked gas which flows 
away from the shock be Vd. Then the distance covered by 
shocked gas passes in a time t is L = Vdt. Therefore, if the iron 
line flux ratio R is known from observations at distance L from 
the shock front then using the function R(t) we can derive the 
value of the downstream velocity Vd = L/t(R), where t(R) is the 
inverse function for R(t). 

The Mach number of the shock and the downstream velocity 
are related by (e.g. [Landau & Lifshitzl l959) 



(2 + ( y -l)M 2 
2yM, 2 -(y-l) 



(22) 



where = Vd/c2 and C2 is the velocity of sound in the post- 
shocked gas. 

On the observational side, it will be important to derive the 
flux ratio of the FeXXV and FeXXVI iron lines from the region 
between the shock front and the considered distance L which 



Fig. 6. Dependence of the Mach number of the shock M on the 
iron line flux ratio R. 



shocked gas covers in a time t = L/V<x, Since non-equilibrium 
ionization can occur at linear scale < 50 kpc (see Sect. 2), we 
choose L = 25 kpc. The flux ratio of the iron lines FeXXV and 
FeXXVI from this region is then 



R 



Jo " (»+24g FeXXV + »+25Q'R R 2 FeXXV ) x dt 

So^ (" + 2 5 G F ; XXV I + "+264r 2 FeXXVl) X dt 



(23) 



Using Eqs. d22l) and d23l . we find the Mach number M of the 
shock as a function of the iron line flux ratio R. The dependence 
M(R) is plotted in Fig. [6] 

Thus, we conclude that the Mach number of the shock can 
be derived from the iron line flux ratio. Methods based on the 
Rankine-Hugoniot jump conditions and on measuring the iron 
line flux ratio are independent for deriving the Mach number of 
the shock. 



6. Conclusions 

The currently operating X-ray imaging observatories provide us 
with a detailed view of the intracluster medium in galaxy clus- 
ters. Chandra's 1" angular resolution, the best among the cur- 
rent X-ray observatories, corresponds to linear scales < 1 kpc at 
z < 0.05 and « 4 kpc at z = 0.3 (the redshift of the 1E065756 
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cluster). This enables us to study hydrodynamic phenomena in 
galaxy clusters, such as shock waves and hot plasma bubbles. 

Metal observations are always limited by the number of X- 
ray photons. For diffuse low surface brightness objects, like 
galaxy clusters, the effective area is a major issue, therefore 
for a reasonable observation time, Chandra metal abundance 
maps will have lower spatial resolution than XMM-Newton (see 
IWerneretal.1 2008). IXO0is planned to be a follow -up mission 
of XMM-Newton and will have a sensitivity much higher than 
XMM-Newton. The expected effective area of the IXO mirror 
and focal plane instruments showing the large improvement at 
all energy (including the 6-7 ke V band) in comparison with those 
of current X-ray observatories is plotted in a figure 0. The larger 
effective area in the 6-7 keV band the higher accuracy of the iron 
line flux measurements is achieved. With IXO we will be able to 
resolve for the first time the metal distribution in the ICM on the 
scales of single galaxies in nearby clusters (simulations of metal- 
licity maps which will be provided by the next generation X-ray 
telescope are given by |Kapferer et al.| 2006). Therefore, the sen- 
sitivity of IXO will provide metal observations near shock fronts 
and in hot plasma bubbles. 

We have considered in this paper the non-equilibrium ioniza- 
tion at linear scales < 50 kpc in galaxy clusters. The necessary 
condition on the baryonic overdensity (see Eq. [8} for the exis- 
tence of non-equilibrium ionization in regions of galaxy clusters 
where 5 > 200 holds for ions of iron. This is because the iron 
atomic number Z=26 is high enough for the threshold value of 
the overdensity, which is proportional to Z 3 (see Eq. [8}, to be- 
come higher than the mean cluster baryonic overdensity. 

The dependence of the He-like and H-like ionic fractions of 
iron on time is given in Sect. 3. We found that the He-like ionic 
fraction of iron does not achieve its equilibrium value during 
the age of the hydrodynamical structures and non-equilibrium 
ionization takes place. 

We calculated the reduced emissivity in the He-like iron 
spectral lines and concluded that the iron abundance derived 
from the assumption that ionization states are in equilibrium pre- 
dicted to be higher than the correct iron abundance value (see 
Fig.©. 

We found that the slowest relaxation process corresponds to 
the increase in the FeXXVII ionic fraction due to decreases in 
the FeXXV and FeXXVI ionic fractions. However, the decrease 
in the FeXXV ionic fraction is much higher than the decrease in 
the FeXXVI ionic fraction during this relaxation process. 

A new way to derive the Mach number of a shock based on 
measurements of the flux ratio of the FeXXV and FeXXVI iron 
lines is proposed in Sect. 5. The advantage of this method with 
respect to the method based on the Rankine-Hugoniot jump con- 
ditions is that the first is more accurate. Fortunately the iron line 
flux ratio is constrained without the effect of hydrogen column 
density (Nh) uncertainties. In practice, the X-ray data can be fit- 
ted in a narrow band containing the FeXXV and FeXXVI lines, 
where the absorption is negligible (see INevalainen et al.1 2009). 
The drawback is that the number of photons is small in this nar- 
row energy band, but the next generation X-ray telescope IXO 
with larger effective area overcomes this drawback and will be 
able to measure the flux ratio of the iron K a lines and, therefore, 
the Mach number of a shock with high precision. Using the nar- 
row energy band instead of the full X-ray spectrum minimizes 
the dependence on calibration accuracy (see INevalainen et al.l 
2003), therefore the FeXXV to FeXXVI lines are insensitive to 



the details of the effective area function in contrast to the con- 
tinuum spectrum. Note that the method based on the Rankine- 
Hugoniot jump conditions uses the densities and temperatures 
derived from the continuum spectrum. 

Another advantage of the proposed method is that it permits 
us to determine independently the Mach number of a shock by 
using measurements of the iron line flux ratio at different dis- 
tances from a shock (see Sect. 5) since it takes into account an 
evolution of ionization states. 

The effect of the apparent iron overabundance under the 
assumption of ionization equilibrium and the slow process of 
helium-like ionic fraction relaxation should be analyzed in 
galaxy clusters by means of future X-ray observatories and 
may have implications in different astrophysical plasmas (e.g. in 
supernova remnants). New high-spectral-resolution instruments 
with higher sensitivity, such as IXO, are needed to measure the 
flux ratio of the iron Ka lines with the purpose of independent 
by determining the shock parameters. 
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